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A new, generalized mean-square algorithm is presented to adjust the 
taps of a transversal adaptive equalizer. Any knowledge of the channel or 
the signaling format can be taken into account and will speed up the con- 
vergence process. The main applications are seen in partial-response sig- 
naling, where the new algorithm eliminates the interaction between the 
individual tap increments. This is achieved by decorrelating the compo- 
nents of the gradient in a fixed weighting matrix prior to the adjustments. 
Convergence is then extremely fast. If the channel has only phase distor- 
tion, a single iteration is sufficient to obtain the optimum tap vector (for 
any timing and carrier phase). This is verified by computer simulations. 
Finally, the new equalizer is compared with another structure recently 
proposed, and some advantages of the new system with regard to imple- 
mentation and flexibility are pointed out. 

I. INTRODUCTION 

Together with sophisticated modulation techniques and bandwidth- 
conserving signal designs, automatic transversal equalizers are the 
prime characteristics of the new generation of high-speed data modems 
that have appeared on the market in the past few years. These modems 
allow data rates in the range of 3600 to 9600 b/s to be transmitted 
over voice-band channels. In addition to point-to-point transmission 
between computers and data processing centers, a large demand for 
high-speed data sets exists in multiparty polling systems. Because a 
message in such a system may consist of only a few hundred bits, it is 
essential that the start-up time of the modem be very short. Prefer- 
ably, the modem start-up time should waste fewer bits than are 
contained in an average message to guarantee a reasonable system 
throughput. This would basically imply that start-up time should be 
inversely proportional to the data rate; this is a very contradictory 
requirement, since more accurate and complex operations are usually 
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required during start-up of a high-speed modem than with a low-speed 
version. Timing recovery, carrier recovery (if required), and training 
of the automatic equalizer are the most important start-up operations. 
The time that must be allowed for the equalizer training is usually 
the major delay in start-up; thus, it seems worthwhile to concentrate 
some of our efforts in this direction. 

During the past few years, much attention has been given to partial 
response signaling. 1-6 Some of these waveforms are particularly 
attractive for data transmission over band-limited channels, since they 
have a spectral zero at the Nyquist frequency. These signals can easily 
be generated with virtually zero excess bandwidth, and yet their time 
response decays faster than that of a sin (x)/x pulse. These advantages 
result from a controlled amount of intersymbol interference in these 
signals, which increases the number of levels in the baseband eye 
pattern. In some formats, the impulse response has odd symmetry 
and the spectrum is zero at dc (which is of particular interest if ssb 
modulation is used). 5 If a traditional transversal automatic equalizer 
with an ms (mean-square) gradient algorithm 67 is used for partial 
response signals, the initial convergence will be very slow. Chang has 
shown that the spread of the eigenvalues of the signal autocorrelation 
matrix can become very large for Class IV partial response signals, 
and he proposed a new equalizer structure that eliminates this problem 
and thus can give fast convergence. 8 Of course, with the appropriate 
modification of some parameters, his conclusions can be extended to 
other partial-response signals as well. 

In this paper, we propose an alternate solution to the problem of 
fast partial-response signal equalization. Instead of changing the 
equalizer structure, we propose to modify the tap updating algorithm. 
The coefficients associated with the new algorithm are given by a 
simple matrix inversion. 

The advantages of the new equalizer structure are not limited to 
partial-response signaling. Any knowledge of the channel or the 
signaling format can be used to speed up the convergence process. 
The classic transversal equalizer evolves as a special case (Nyquist 
signaling) of the more general configuration that is proposed in this 
paper. 

II. A GENERAL MEAN-SQUARE ALGORITHM 

Let us assume a transversal equalizer structure with N taps, with a 
tap signal input vector x* and a tap coefficient vector c at time t Q + kT. 
The equalizer output y* is given by 

yu = tfc (1) 
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and will not, in general, be exactly equal to the desired reference value 
dk, the error being 

e k = Vk — d k . ( 2 ) 

The ms error that appears at the output will thus be 

S = E{el\, (3) 

and we now choose a tap-updating algorithm that adjusts c in such 
a way that the ms error (3) is minimized. The gradient of e 2 with 
respect to c can easily be obtained by combining (1) to (3), 

g = - 2 % = E{ekX k } =Ac-v, (4) 

where we have introduced the signal autocorrelation matrix, 

A=E[x*Z), ( 5 ) 

and the signal correlation vector, 

v = E{x k dk\. ( 6 ) 

The optimum tap vector in the mean-square sense that minimizes (3) 

is obtained when we set the gradient to zero, 

Copt = A~ 1 y. (7) 

We propose a tap-updating algorithm of the form 

C m+ 1 = C m - Qmgm, (8) 

where g m is the gradient evaluated with tap-vector c m after the wth 
iteration, and Q m is a nonsingular matrix, which is discussed below. 
Obviously, an algorithm of the general form (8) will stop updating 
(turn itself off) when c = c opt . If the algorithm converges at all, it 
will thus converge to the ms solution (7). The traditional steepest- 
descent gradient algorithm is obtained as a special case from (8) 
if we set Q m = 01 (or a more general diagonal matrix instead of the 
unity matrix / if tapering should be included). 
The algorithm (8) can be written in the form 

c m+ i = (/ - Q m A)c m + QmV. (9) 

If we introduce the tap error vector 

ifm = Cm C pt> 

it is easy to show that 



(10) 



hl -(I- Q m A)i, m = n (/ - Q«A)+i. (11) 



=1 
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Obviously, the convergence can be controlled to some degree by the 
selection of Q m , since this choice affects the eigenvalues of / — Q m A. 
This is discussed in Section III. 

III. SELECTION OF THE Q MATRIX 

First we will make some general remarks without specifying that 
the signal format be partial response. If the channel characteristics, de- 
modulating carrier phase, and timing instant were perfectly known, 
we would select Q = Q m = A~ l and then obtain c opt after the first 
iteration. Actually, if that knowledge were available, it would, of 
course, be easier to preset the taps to c op t or to introduce a suitable 
fixed compensation network. Unfortunately, the transmission channel 
is not completely known. If, however, the spread within these channels 
is not extreme, we can reduce the average equalizer training time if we 
select a good estimate of the optimum tap vector for initial presetting. 
For example, 

ci = EUr^t], (12) 

where the expected value is taken over the ensemble of channels, and 
a similar estimate for Q, 

Q m = (}JjE{Ai}y-\ (13) 

For a particular channel with 

Ai = E{A i \ - A., (14) 

the critical matrix Bi = I — Q m Ai then becomes 

B { = (1 - ft.)/ + m lElAi}J-*Ai. (15) 

If we select /3 m = 1, the average Bi matrix is equal to the zero matrix. 
For a randomly selected channel, the B t matrix will, on the average, 
be much closer to the zero matrix than with a gradient algorithm, and 
faster convergence may be obtained. Modern digital signal-processing 
techniques allow quick and easy change (or selection) of (12) and (13) 
with different read-only memories at a fraction of the costs that are 
involved when traditional compromise equalizers have to be changed- 
Furthermore, even different signal formats may be processed effi- 
ciently with the same modem. 

IV. APPLICATIONS TO PARTIAL-RESPONSE SIGNALING 

Our results have so far been quite general, as we have not yet 
defined a special signaling technique. As was already mentioned in 
Section I, we see a particular promising application to processing zero 
excess-bandwidth partial-response signals. For such signals, the A 
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matrix depends only on the amplitude characteristics of the channel 
and is independent of the phase parameters 8 (timing phase, carrier 
phase, and phase characteristic of the channel). Since the amplitude 
response of voice-grade channels is relatively constant over the trans- 
mission band (or at least under much better control than the phase 
response), we may select an estimate 

Q m = ftJLo 1 , (16) 

where A is the ideal partial-response signal matrix. More specifically, 
for Class IV signaling, the matrix elements of A are 

[Ao]ik = S« ~~ 2^i,t+2 — h^i.t-2 (1«) 

and for duobinary are 

[Ao]« = «« + !«<.*+! + !*.\*-i- ( 18 ) 

The inverse of both (nonsingular) matrices may be expressed as 

Ao 1 = (I-G)- 1 = I + G + G 2 + • • -, (19) 

which shows that the elements of A _1 are rational numbers, and are 
symmetric with respect to both matrix diagonals. 

Having selected an ideal channel as an estimate for Q, our algorithm 
takes the form 

♦«+i = (J - & m Ao l A)if m = B m $ m . (20) 

Obviously, if we have only phase distortion, then A = A and we 
may select /3 m = 1. The equalizer will then reach its optimum tap 
vector within a single iteration because the matrix B n is zero. 

If the channel has amplitude distortion, then A ?± A and thus 
B m ?± 0. The eigenvalues X, and eigenvectors w, of A^ l A are defined by 

A^Awt = X,w,-. (21) 

Following a procedure shown by Chang, 8 we premultiply both sides by 
vr\A o and obtain 

Xj = ^i. (22) 

Since both A and A are symmetric, positive, definite matrices, we 
conclude that all X,- are positive. Chang 8 has further shown that the 
eigenvalues can be bounded by 

|ff(-)|i*£X«3! I* MIL* (23) 

where H(u) is the quotient between the actual and the estimated 
channel transfer function. For convergence, the spectral radius p of 
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B m must be less than unity ; this is guaranteed if we choose 

< p m < ^- , (24) 

since the eigenvalues m of B m are given by 

Ml = 1 - pnki. (25) 

If the amplitude distortion is small, A^A is close to the identity 
matrix, and the eigenvalues X, will be closely scattered around unity. 
We would then select /3 m tt const tt 1, and the matrix B would be 
close to the zero matrix with all eigenvalues m close to zero. The tap 
error vector after m iterations is given, from (20), 

4? m = (I - PA^A)^ = B m $ (26) 

if /3 is constant (this is, of course, replaced by the usual matrix product 
if P = /3 m ). The squared magnitude of § and the excess MS error are 
given by the expressions 

|<M 2 = <tf£ 2 "**o (27) 

e 2 = $B m AB m $ , (28) 

and can be used to estimate the convergence speed. The procedure 
and the results are practically identical to those of Chang 8 and are, 
therefore, not repeated here. 

V. IMPLEMENTATION 

Figure 1 is a basic block diagram of the new equalizer. A version with 
only three taps is shown for clarity, although a much larger and prefer- 
ably even number of taps would actually be used. The upper part 
represents a traditional ms equalizer with the well-known circuitry 
for correlating the error signal with each tap signal to obtain the 
gradient g = Ac — v. Because of the particular structure of A for 
partial response signals, the components of g are highly correlated. 
The tap corrections obtained with a direct gradient algorithm would 
thus not be independent of each other. This interaction (which would 
cause a very slow convergence) can be eliminated if the gradient is 
"decorrelated" in the fixed weighting matrix Q shown in the lower 
part of Fig. 1. The resulting outputs are now decoupled and can be 
used to adjust the taps free of interaction to obtain rapid convergence. 

A somewhat different arrangement is shown in Fig. 2. Here the 
(same) Q matrix is placed between the tap signals and the correlators. 
It is obvious that this arrangement is equivalent to Fig. 1 because the 
error signal is common to all correlator inputs. 

148 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1975 



SIGNAL 
INPUT 



0. MATRI) 



- * * i* 

g g g 



L® Up La. 



OUTPUT 



■REFERENCE 



~--. CORRELATORS 
(MULTIPLY + AVERAGE) 



Ac 2 



Ac 3 



CONTROL SIGNALS 
(-FOR TAP 
INCREMENTS 



Fig. 1 — Use of Q matrix in a first implementation of the new algorithm. 
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Fig. 2 — Use of Q matrix in a second implementation of the new algorithm. 
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Fig. 3 — Use of P matrix to obtain a modified equalizer structure. 

As a comparison, we show in Fig. 3 a three-tap equalizer proposed 
by Chang. 8 Here the input signal is first transformed by a weighting 
matrix P into a set of orthonormal components that are then processed 
as in a conventional equalizer. It is, in fact, identical to a conventional 
equalizer except that the input signal vector x is replaced by Px, 
where P is a suitably chosen matrix such that P T P = Aq 1 . Initial 
presetting of the equalizer and constraining the range of the variable 
gain coefficients can, however, be quite complicated with this structure. 
These difficulties are avoided in the new proposal, since the equalizer 
structure remains unchanged and only the updating algorithm is 
modified. 

Although the number of multiplications is the same in all three 
structures (Figs. 1 to 3), we would like to point out one very important 
difference: in Chang's structure, the accuracy of the transformation 
is very critical. The input signals must be represented by their full 
precision (say, 10 bits). To avoid round-off errors when all the products 
are summed up, more than 12 accurate bits of these products must be 
processed and can then (after summation) probably be rounded off 
to 10 bits for the transformed tap signal. Any errors introduced in this 
process will result in undesirable output noise that cannot be compen- 
sated for by the equalizer. In the new equalizer, however, only the tap 
increments are affected by errors in the matrix calculation. Obviously, 
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a much lower accuracy is needed; in the extreme case, the correct 
polarity will be sufficient. Even with comparatively large errors, the 
tap-vector noise can be kept low by suitable scaling of the increments. 
We thus conclude that the new equalizer has inherent advantages as 
far as complexity, cost, size, speed requirements, and time-sharing 
possibilities are concerned. Furthermore, simpler algorithms can be 
derived from this one. If, for example, the structure of Fig. 2 is used 
for an algorithm with tap signals quantized to one bit, the matrix 
operation is reduced to controlling the addition and subtraction of 
the coefficients and may, for a medium-sized equalizer, be cheap and 
efficiently replaced by a suitably programmed read-only memory. 
For such nonmean-square algorithms, the optimum Q matrix will, in 
general, of course, be different from Aq 1 . 

A few words should be said about the structure of the matrix 
Q = A _1 in the special case of partial-response Class IV signaling. It 
can be shown that every second diagonal array is zero ; furthermore, 
if N is even, every second row (or column) is a shifted version of the 
previous one. These properties are identical to those of Chang's P 
matrix (the remaining nonzero coefficients are, however, different from 
those of the P matrix) and can be used to simplify the signal-processing 
operations. 

Finally, we would like to point out that an all-digital implementation 
for higher data speeds (group band) would still be extremely complex 
and expensive. For such applications, an array of operational amplifiers 
with suitable weighting resistors (probably mounted on a ceramic) 
would be a far more economic solution with the present state of the art. 



VI. SIMULATIONS 

To test the proposals made in this paper, a 15-tap ms equalizer was 
simulated. A 4800-b/s, ssB-modulated, partial-response, Class IV 
signal was selected. A parabolic delay characteristic was assumed 
with 1-ms delay difference between the passband channel center and 
edges. Table I shows the coefficients of the matrix Q = A? 1 . In Fig. 4, 
the equalizer convergence is shown for four different timing phases 
that are separated by 90 degrees each. After each iteration, the ms dis- 
tortion of the overall impulse response is plotted. It is seen that the 
minimum is basically achieved with the first iteration (the small resid- 
ual error after the first iteration results from the time-domain trunca- 
tion and round-off errors in the components of the initial impulse 
response during the simulation). Some additional simulations were 
made with an ms stochastic approximation algorithm (adjustments at 
the symbol rate) and with amplitude distortion; they also showed a 
clearly improved convergence when the decorrelation matrix was used. 
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Table I — Decorrelation matrix for N = 15 
(elements multiplied with 36) 
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Fig. 4: — Convergence of new algorithm for different timing phases. 
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VII. CONCLUSIONS 

The spread of the eigenvalues of the signal autocorrelation matrix 
can become very large for signals with nonflat amplitude spectra. If 
the gradient with respect to the tap gains of an ms equalizer is calcu- 
lated, its components are highly correlated. The resulting interaction 
of the individual gain adjustments may cause a very slow equalizer 
convergence. We have proposed a new ms algorithm that eliminates 
these problems. With this new algorithm, a set of uncorrelated tap 
increments is calculated so that the individual adjustments are de- 
coupled from each other. This can be achieved in either of two ways : 
The gradient can first be calculated in the conventional way (corre- 
lation of error signal and tap signals) and then be multiplied by a 
decorrelation matrix Q to produce the final correction vector. On the 
other hand, tap signals can first be passed through such a matrix and 
the resulting outputs can then be correlated with the error signal. The 
matrix Q is chosen as the inverse of the expected (average) signal 
autocorrelation matrix A. The new algorithm was simulated for partial 
response signals with zero excess bandwidth. In the case of a flat 
channel, the optimum tap coefficients can then be determined within 
a single iteration, independent of delay distortion, carrier phase, and 
timing instant. This was confirmed by simulation. Compared with 
other recent proposals for fast start-up, the new algorithm features 
simpler presetting and requires less accuracy in the matrix trans- 
formation. Furthermore, this method can easily be extended to systems 
with nonmean-square cost functions. 
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